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Recent deep Chandra LETG+HRC-S observations suggest the possibility that RX 
J1856. 5-3754 is a compact star whose radiation radius is 3.8-8.2 km. In this paper, we 
systematically calculate the mass-radius relations of quark stars within the bag model. As- 
suming that RX J1856. 5-3754 is a pure quark star, we derive an upper limit on its mass for 
various sets of bag-model parameters. 
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§1. Introduction 

The possibility that compact stars supported by the degenerate pressure of quark 
matter exist has been investigated by many authors (see, e.g., Refs. l)-20)). The 
Qh' most commonly accepted picture of such quark stars is that a star containing quark 

matter in its core region, surrounded by hadronic matter, might exist in the branch of 
neutron stars; 13 ) this hypothesized type of star is often referred to as a "hybrid star" . 
If hadronic matter in neutron stars can undergo a strong first-order phase transition, 
it is possible that a third family of more compact stars exists. 21 ) Such a third family, 
. ^ | arising from a deconfinement phase transition, is predicted by a very restricted class 

^ ■ of models of quark and hadronic matter. 14 )' 15 ) As suggested by Witten, 10 ) there is 

another possible form of quark stars: If the true ground state of hadrons is "strange 
matter" , bulk quark matter consisting of approximately equal numbers of u, d, and 
s quarks, then self-bound quark stars (or "strange stars") could occur at masses and 
radii on the order of or even smaller than those of typical neutron star, ~ 10 km and 
~ IAMq. In the absence of reliable information about the equilibrium properties of 
hadronic and quark matter at high densities, however, it is impossible to tell which 
kind of quark stars is favored. 

Recently, Drake et al. 22 ) reported that the deep Chandra LETG+HRC-S obser- 
vations of the soft X-ray source RX J1856. 5-3754 reveal an X-ray spectrum quite 
close to that of a blackbody of temperature T = 61.2 ± 1.0 eV, and the data contain 
evidence for neither spectral features nor pulsation. *) They argued that the derived 



*' They placed a 99% confidence upper limit of 2.7% on the unaccelerated pulse fraction from 
10 -4 to 100 Hz. Burwitz et al. 25 ^ reported that a subsequent observation of RX J1856. 5-3754 with 
XMM-Newton allows the upper limit on the periodic variation in the X-ray region to be reduced to 
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interstellar medium neutral hydrogen column density is Njj = (0.8-1.1) x 10 20 cm~ 2 , 
which, together with the results of recent HST parallax analyses, yields an estimate 
of 111-170 pc for the distance D to RX J1856. 5-3754. Combining this range of D 
with the blackbody fit leads to a radiation radius of Roq = 3.8-8.2 km, which is 
smaller than typical neutron star radii and thus suggests that the X-ray source may 
be a quark star. 

In a subsequent work, Walter and Lattimer 23 ) showed that the blackbody model 
of Drake et al. does not reproduce the observed UV/optical spectrum. 24 ) They 
succeeded in fitting the two-temperature blackbody and heavy-element atmosphere 
models developed in Ref. 24) to the observed spectrum, ranging from X-ray to op- 
tical wavelengths. They found that the radiation radius inferred from such a fitting 
lies between 12 km and 26 km, which is consistent with that of a neutron star. How- 
ever, this model is not effective for explaining the lack of spectral features. Most 
recently, Braje and Romani 26 ) suggested that a two-temperature blackbody model, 
which can reproduce both the X-ray and optical-UV spectral data but seems to be 
inconsistent with the fact that pulsation is not detected, is indeed compatible with 
the interpretation that the object is a young normal pulsar, whose nonthermal radio 
beam misses Earth's line of sight. However, this model cannot answer why there are 
no features in the observed X-ray spectrum. In the absence of a model that accounts 
for all the observational facts, in this paper we employ the simple picture of Drake 
et al. 22 ) based on the uniform temperature blackbody fit to the X-ray spectrum. 

From the radiation radius of RX J1856. 5-3754 inferred from the blackbody fit to 
the X-ray spectrum, the upper bound of the star's mass can be derived. 27 ) -29 ) This 
is because the radiation radius is larger than the true radius by a factor arising from 
gravitational redshift effects. The analysis of Drake et al. yields an upper bound of 
0.5-1.1Mq. In previous studies, 27 )' 29 ) the mass-radius relations of pure quark stars 
were calculated from several kinds of equations of state (EOS) for quark matter, 
and these results were compared with the constraints on the mass and radius of RX 
J1856. 5-3754 derived from the inferred radiation radius. In this framework, details 
of the relation between the mass of RX J1856. 5-3754 and the EOS remains to be 
clarified. Further systematic investigations allowing for uncertainties in the EOS arc 
needed; it is expected that these uncertainties can be reduced using Monte-Carlo 
calculations of lattice gauge theory at non-zero density, though such calculations at 
present are beset by technical problems. 

Taking into account the above points in this paper, we consider the EOS of quark 
matter in the space of the parameters characterizing the bag model and investigate 
the question of how massive a quark star could be, given the inferred radiation radius. 
In §2, we construct quark stars using the bag-model EOS of quark matter. Section 
3 is devoted to derivation of the mass-radius relation and its application to the soft 
X-ray source RX J1856. 5-3754. Conclusions are given in §4. We use natural units 
in which h = c = &b = 1 throughout the paper. 



1.3% at 99% C.L. from 10" 3 to 50 Hz. 
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§2. Quark star models 



To construct stars composed of zero-temperature uds quark matter, we start 
with the thermodynamic potentials, f2 q , for a homogeneous gas of q quarks (q = 
u,d,s) of rest mass m q to first order in a c , where a c = gl/^n is the fine structure 
constant associated with the QCD coupling constant g c . In Refs. 30), 31) and 9), 
expressions for Q q are given as sums of the kinetic term and the one-gluon-exchange 
term at the renormalization scale A = m q . We assume that m u = = and write 
f2n as 



with 
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F(x s ) = x s 7] s - In (x s + f] s ), 



(2-4) 



where x s = fj, s /m s , rj s = \Jx'j — 1, (i q is the chemical potential of q quarks, a c = 
a c (A), and m s = m s (A). In this paper we choose A = fi s .*^ 

In a star, electrons are distributed in such a manner that each local region is elec- 
trically neutral. (We ignore muons, which are more massive than electrons.) Since 
generally the electrons are dense enough to behave as an ultrarelativistic ideal gas, 
we can write the electron thermodynamic potential Q e in the massless noninteracting 
form, 



12vr 2 ' 



(2-5) 



where [i e is the electron chemical potential. 

In the bag model, we express the total energy density p as 32 ) 

i=u,d,s,e 



(2-6) 



where B is the bag constant, i.e., the excess energy density effectively representing 
the nonperturbative color confining interactions, and rii is the number density of i 



*•* Farhi and Jaffe 32 ' chose the renormalization scale as A = 313 MeV in describing self-bound 
uds quark matter in bulk; in such matter, the one-gluon-exchange interaction energy is small com- 
pared with the kinetic energy. Here we describe uds quark matter forming a system of stellar size 
with a much wider density range. The present choice (A = fi s ) is sufficient to keep the interaction 
energy relatively small for such a density range. 
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particles, given by 

n ■ — — 

dfii ' 

Then, we obtain the pressure p as 



Ui = (2-7) 



p = HBUB ~ P 

= - ]T Qi-B, (2-8) 

i=u,d,s,e 

where ps = J2i= u ,d,s,e n if J 'i/ n B is the baryon chemical potential, and 

nB = g ( n u + n d + n s ) (2-9) 

is the baryon density. 

The values of the parameters m s , a c and B were obtained from fits to light 
hadron spectra in Refs. 33) -35) (see Table I). However, the values yielded by such 
fits do not correspond to bulk quark matter. 32 ) Allowing for possible uncertainties, 
we consider values in the ranges 130 < B 1 / 4 < 250 MeV, < m s < 300 MeV, and 
< a c < 0.9. The ranges of m s and a c so chosen are not consistent in all cases with 
the values in Table I, but they are consistent with the Particle Data Group 36 ) values 
renormalized at an energy scale of interest here. 



Table I. Bag-model parameters fitted to hadron mass spectra. 



B 1/4 (MeV) 


nis (MeV) 


Ctc 


Reference 


145 


279 


2.2 


T. DeGrand et al. (1975) 33) 


200-220 


288 


0.8-0.9 


C.E. Carlson et al. (1983) 34) 


149 


283 


2.0 


J. Bartelski et al. (1984) 35) 



To obtain the equilibrium composition of quark matter in its ground state at 
a given baryon density or chemical potential, the conditions of equilibrium with 
respect to the weak interaction and of overall charge neutrality must be imposed. 
These conditions are expressed by 

Pu = Pd-^e, (2-10) 
Pd = Ps, (2-11) 

and 

2 11 . 

2 n u ~ ~ = n e- (2d2) 

From Eqs. (2-l)-(2-12), we can obtain the energy density and the pressure as func- 
tions of tib (or pb)- 

In Figs. 1 and 2, we plot the energy per baryon (p/ris) as a function of n^ 1 . 
We note that p/ris becomes large as the strange quark mass m s increases, because 
p s is an increasing function of m s for fixed n s . We also find that p/ns increases as 
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Fig. 1. Energy per baryon of electrically neutral quark matter in j3 equilibrium (minus the nucleon 
rest mass mjv) as a function of rig 1 , calculated for a c = and _B 1//4 = 145 MeV. The dotted 
curve is the result for ud quark matter (m s = oo). The dashed (long dashed) curve is the result 
for uds quark matter with m s — 300 (150) MeV. The solid curve represents the result for a 
massless three-flavor system (m s — 0). 



a c increases. This is because the one-gluon-exchange interaction is repulsive, due to 
the contribution from the exchange of transverse (magnetic) gluons in a relativistic 
quark plasma. 

Figure 3 depicts the equation of state of ground-state quark matter, i.e., the 
pressure as a function of the energy density. From this figure, we see that the EOS 
becomes softer as B increases and hence as the color confining force increases. Such 
a feature is explicitly expressed in Eq. (2-8). We also find that the EOS becomes 
softer as s quarks become more massive. This softening stems from the fact that a 
nonzero but small m s acts to reduce the kinetic pressure of s quarks, as can be seen 
from Eq. (2-3). We remark in passing that as p increases the quark matter EOS 
becomes closer and closer to the massless ideal-gas limit (p = p/3). 

We proceed to calculate the structure of nonrotating quark stars as in Refs. 11) 
and 12). Such calculations can be performed by incorporating the EOS models 
obtained above into the general relativistic equation of hydrostatic equilibrium, i.e., 
the Tolman-Oppenheimer-Volkoff (TOV) equation, 37 ) 

dp(r) [p( r ) + p( r )] [M(r) + 47rr 3 p(r)] /() 

dr 

r 



2GM(r) 
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Fig. 2. Same as Fig. 1, but a c = 0.6. 
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. 3. Equation of state of uds quark matter in the ground state, calculated for B 1 ^ 4 = 145, 200 
MeV and a c — 0.6. The dashed, long dashed, and solid curves are the results for m 3 — 300, 150, 
MeV, respectively. The long dashed curve represents the EOS of a massless ideal gas (p = p/3). 
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where G is the gravitational constant, r is the radial coordinate from the center of 
the star, and M(r) is the gravitational mass of the stellar portion inside a surface of 
radius r, which can be obtained by integrating the equation for mass conservation, 

^=Wp(r), (2-14) 
dr 

from to r. Then we can determine the radius, R, of the star from the condition 
that the pressure becomes zero, i.e., p(r = R) = 0. For various values of the energy 
density, po, at the center of the star, we finally obtain the radius R and the mass 
M = M(R). 
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Fig. 4. Masses of pure quark stars with central energy density po, calculated for B 1 / 4 = 145,200 
MeV and a c = 0.6. The dashed, long dashed, and solid curves are the results for m 3 — 300, 150, 
MeV, respectively. 

In Fig. 4 we plot the mass M as a function of p$. In this figure, we include the 
M-po relation up to the maximally allowed value of the mass; at this point, the star 
is no longer gravitationally stable. 38 ^ We also exhibit the radius R as a function of 
po in Fig. 5, where the right end of the R-po relation corresponds to the maximum 
mass for a quark star As can be seen in both Figs. 4 and 5, the radius and the 
mass become smaller for fixed po as B and/or m s increases, which leads to softening 
of the EOS (see Fig. 3). In other words, an increase in B and/or m s for a fixed 
radius or mass tends to increase pq. 

For m s — 0, the mass, radius and central energy density of a stable quark star of maximum 
mass are M = 2.00M© {B /B) 1/2 , R = 11.1 km (B /B) 1/2 and p = 1.91 x 10 15 g cm" 3 (B/B ). 10) 
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Fig. 5. Same as Fig. 4, except that the quantity plotted is the radius rather than the mass. 



§3. Upper limit on the mass of the possible quark star 

In this section we derive the mass-radius relations of pure quark stars for various 
bag-model parameters and compare these relations with the inferred radiation radius 
of the soft X-ray source RX J1856. 5-3754. We can straightforwardly calculate the 
mass-radius relations from the M-po and R-po relations obtained in the previous 
section. The results for B 1 ^ = 145, 200 MeV, m s = 0, 150, 300 MeV, and a c = 0.6 
are plotted in Figs. 6 and 7. The M-R relation ends at the maximum mass, where 
the gravitational instability sets in, while the mass M, when small, behaves as ~ R 3 , 
due to the vacuum pressure B. 

We now compare the obtained mass-radius relations of quark stars with the 
inferred radiation radius of RX J1856. 5-3754. According to the analysis of the 
deep Chandra LETG+HRC-S observation of RX J1856. 5-3754 carried out by Drake 
et al., 22 ) the X-ray spectrum is quite close to that of a blackbody of temperature 
T = 61.2±1.0 eV, and the X-ray luminosity is L x ~ 6 x 10 31 (L>/140 pc) 2 ergs" 1 . 22 ) 
Then, we obtain the relation, Lx = ^R^ (7r 2 /60T 4 ) between the luminosity and 
the temperature with the radiation radius Roo given by 

R 

Roo = j (3-1) 
' 2GM y ' 

R~ 



The factor — 2GM/R represents the redshift effect in a strong gravitational 
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Fig. 6. Mass-radius relations of pure quark stars, calculated for a c = 0.6 and _B 1//4 = 145 MeV. 
The dashed, long dashed, and solid curves are the results for m s — 300, 150, MeV, respectively. 
The light shadowed region is allowed by the inferred radiation radius, R x — 3.8-8.2 km, of RX 
J1856. 5-3754. The dark shadowed region is excluded by the condition that R is larger than 
2GM. 
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Fig. 7. Same as Fig. 6, but for B 1/4 = 200 MeV. 
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field. Drake et al. 22 ) inferred i?oo = 3.8-8.2 km from the distance to RX J1856.5- 
3754, D= 111-170 pc. 

Given such an inferred radiation radius Roo, we can obtain information about 
the radius R and the mass M of the star from Eq. (3-1). In Figs. 6 and 7, the light 
shadowed region denotes the region of R and M allowed by the radiation radius 
Roo = 3.8-8.2 km. The dark shadowed region in these figures is excluded by the 
condition that the radius R of a radiation emitter be larger than a black hole surface 
of radius 2GM. From these figures, we can see that the mass upper limit, M up , 
allowed by R^ increases from ~ OAMq to ~ 1-Mq when B 1 ^ increases from 145 
MeV to 200 MeV, while being almost independent of m s . We remark that the mass 
lower limit allowed by R^ is of order O.1M . 

In Fig. 8 we plot the contours of the mass upper limit M up in the space of the 
parameters B 1 ' 4 and m s , calculated for a c = 0,0.3,0.6,0.9. From this figure, we 
find that the contours are peaked around B 1 / 4 ~ 200 MeV. The left side of this 
peak results from the boundary of the allowed M-R region by the inferred radiation 
radius of RX J1856. 5-3754, while the right side is due to the maximum mass of 
gravitationally stable quark stars. We remark that for the smallest values of B, in 
addition to strange matter, even ud quark matter at zero pressure can be more stable 
than normal nuclei, as shown by Farhi and Jaffe. 32 ) We also note that the absolute 
stability of strange matter at zero pressure, which should allow for the existence of 
strange stars to exist, limits the bag-parameter region to the left side of the dot- 
dashed curves. It might be possible to relax this limitation if we consider that the 
bag constant depends on the pressure, as discussed in the Appendix. 

§4. Conclusions 

In this paper we have calculated the mass-radius relations of pure quark stars 
for various sets of bag-model parameter values. Assuming that RX J1856. 5-3754 
is a pure quark star, we have derived an upper limit on its mass, which is of order 
O.5-1M . In such a star of M ~ 1M , the central energy density would be ~ 10 16 
g cm~ 3 . The present systematic analysis essentially includes the cases treated in 
previous analyses utilizing the EOS models based on perturbation theory 16 )' 17 ) and 
finite-temperature lattice data. 18 ) Because this paper treats pure quark stars alone, 
however, many questions remain. We now discuss some of these. 

Unless uds quark matter forms a more stable self-bound system than 56 Fe, a pure 
quark star would likely be unstable with respect to hadronization. The resultant star, 
no matter whether it be a hybrid star or a member of the third family mentioned in 
the Introduction, should have a larger radius for fixed central energy density with 
a hadronic shell surrounding a quark matter core. The extent to which the radius 
would increase depends strongly on the model of dense hadronic and quark matter 
used. Due to this increase, it is possible that the situation in which the bag-model 
parameter sets do not yield the absolute stability of uds quark matter would be in 
contradiction with the inferred radiation radius of RX J1856. 5-3754, except for the 
case of possible parameter sets including large B for which the mass-radius relation 
of the third family enters the region allowed by the inferred radiation radius. 
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Fig. 8. Contours of the mass upper limit, M up , allowed by the inferred radiation radius of the soft 
X-ray source RX J1856. 5-3754 in B 1/4 -m s plane, calculated for (a) a c = 0, (b) 0.3, (c) 0.6, 
and (d) 0.9. In each panel the contours for M up = IMq (solid curves), M up = 0.9Mq (dashed 
curves), M up — 0.7Mq (dotted curves), and M up = 0.5Mq (long dashed curves) are included. 
For reference, we denote by the dot-dashed curves the boundary below which strange matter is 
absolutely stable at zero pressure (for details, see Fig. 1 in Farhi and Jaffe 32 '). 



The possibility for the existence of a quark star of the third family depends 
on the condition that the velocity of sound be drastically enhanced at a deconfine- 
ment transition. 21 ) It may be possible for this condition to be satisfied, as shown in 
Refs. 14) and 15) in the case in which the deconfinement transition is accompanied 
by a quark-hadron mixed phase including uds quark droplets embedded in a gas of 
electrons, in a sea of hadronic matter. This phase may be realizable because the 
presence of s and d quarks acts to reduce the electron Fermi energy and to enhance 
the proton fraction in hadronic matter. However, there are several reasons why it 
might be difficult for such a mixed phase to appear. First, the quark-hadron in- 
terfacial energy, which is poorly known, might be large enough to make the mixed 
phase energetically unfavorable. 39 ^ Second, even for fairly small values of the inter- 
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facial energy, electron screening of quark droplets could act to destabilize the mixed 
phase. 40 ) Finally, even if the mixed phase is energetically favorable, it is uncertain 
how the mixed phase would nucleate in a neutron star. 41 ) 

If uds quark matter is absolutely stable, strange stars could exist, with a pos- 
sible crust of nuclear matter. Such a crust would likely be composed of a Coulomb 
lattice of nuclei immersed in a roughly uniform gas of electrons, but not in a sea 
of neutrons. n ) This nuclear matter including free neutrons would be surrounded by 
absolutely stable uds quark matter; the electric charge of the quark region, which 
is generally positive, would prevent the normal nuclei from contacting the quark 
matter. The effects of the crust on the strange star structure were investigated by 
Zdunik, 19 ) He showed that the maximal presence of the crust would qualitatively 
change the mass-radius relation only for M < O.IMq. Determining to what extent 
the crust would prevail in a strange star requires information about the formation 
and evolution of the star. We remark that even if RX J1856. 5-3754 is a strange star, 
the surface of the star should be hadronic. This is because the observed spectrum 
is of blackbody type, while at the observed X-ray luminosity, the spectrum emit- 
ted from a bare quark matter surface would differ drastically from the blackbody 
spectrum. 42 ) 

Information about the mass of RX J1856. 5-3754 is useful to distinguish between 
the possible forms of compact stars. In order to derive a constraint on the mass from 
the observed X-ray luminosity Lx, it is instructive to note the possible origins of 
Lx ■ Initial cooling and accretion of interstellar material are generally considered to 
be responsible for Lx- 43 ^ We may thus assume that the gravitational energy release 
of the accreting material is smaller than Lx, i.e., Lx ^ GMM /R, with the Bondi 
accretion rate M onto a star moving supersonically in an ambient medium given 
by 44 ) 

(GM\ 2 

Here, A is an O(l) constant related to the EOS of the accreting matter, v is the 
velocity of the star, and pu is the hydrogen mass density (= m^nn with nucleon 
mass rriM and hydrogen number density n#). The velocity v has been determined 
asti~ 200(D/140 pc) km s -1 from observations of the proper motion of the optical 
counterpart of RX J1856. 5-3754. 23 ) In order to estimate nu , we follow the argument 
of Ref. 28). The hydrogen number density is roughly given by njj = Nh/Rh, where 
Nh is the column density to the star derived from the blackbody fit to the X-ray 
spectrum as N H = (0.8-1.1) x 10 20 cm -2 , 22 ' and R H is the size of the local high 
density region. Rh could be on the order of 100 AU or even less, as is known from 
observations of the interstellar medium using the 21 cm Hi line. 45 ) From Eq. (4-1), 
we can thus estimate the mass of the star as 

/ v \ ( D \ 5/3 / N H \ _1 /3 / r h \i/3 
11 " U LUti \ 200 km J iTTolW Uo 20 cm- 2 J [wo AU J 

> \ 1/3 

(4-2) 
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Uncertainties in Eq. (4-2) come mainly from the velocity v and the local high den- 
sity scale Rh- The actual values of v and Rh depend on the poorly known radial 
velocity and density profile in the R CrA molecular cloud believed to contain RX 
J1856. 5-3754. Therefore, M may in fact be much larger than O.4M . However, it 
is interesting to note that the mass range estimated from the X-ray luminosity Lx, 
which is almost independent of the spectral analysis leading to the allowed M-R 
region depicted in Figs. 6 and 7, may play a role in narrowing the allowed M-R 
region and hence in clarifying the composition of this compact star. 

Elucidating the formation of quark stars is another important problem. If quark 
stars of masses smaller than ~ 1M are converted from neutron stars with canonical 
masses of order 1.4M , material in a neutron star should be partly ejected out of 
the star and partly changed into quark matter. Thus, possible scenarios for the 
formation of low-mass quark stars (some of which are discussed in Refs. 46) and 28)) 
depend inevitably on the initial and final configuration of the stars. 
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Appendix A 

Effects of the Pressure Dependent Bag Constant 

In this appendix, we modify the bag model adopted in the main text by allowing 
the bag constant to depend on the pressure. In the model thus modified, it is possi- 
ble for the bag constant inside the star to become effectively larger than that at the 
surface. Then, the absolute stability condition of strange matter at zero pressure, 
which is characterized by the dot-dashed curves in Fig. 8, can be changed signifi- 
cantly. In the main text, we studied the mass-radius relations of pure quark stars, 
with B fixed throughout the interiors, and extrapolated our calculations beyond the 
absolutely stable region. As we see below, however, the stable region can be enlarged 
by the pressure dependence of the bag constant. 

As discussed in Refs. 18) and 20), we can generally parameterize the EOS of 
strange matter in the massless limit (m s — > 0) as 

p=j(p-4B), (A-l) 

where A is a parameter, whose value is three in the standard bag model. Then, we 
can transform Eq. (A-l) into 



P= ^(p-4-Bcff), 



(A-2) 
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where B e s is the effective bag constant, which depends on the pressure p (or the 
energy density p) as 

B cS = ^(A-3)p + B, 




(A-3) 



B e Q reduces to B in the limit p — > 0. From Eq. (A-3) we find that if A is larger than 
three, B e Q becomes larger than B inside the star. In the model used in Ref. 18), A 
takes a typical value of 3.8, while A « 2.2 in Ref. 20). Here we consider the range of 
values A = 2-6. We remark that there should be an upper bound on A; at non-zero 
pressures, if A is far greater than unity, the effective bag constant is increased to 
such an extent that hadronic matter can be energetically favored even at extremely 
high pressures. 

We now consider calculations of mass-radius relations of pure quark stars from 
the EOS (A-l). By repeating the calculations performed in §2 for B e g rather than 
B, we can compute the central pressure of a stable quark star of maximum mass 
for various values of A and B. In Fig. 9, we plot the contours of the effective bag 
constant B l c ^ at this central pressure in the B 1 / i -A plane. We denote by the dot- 
dashed curve the boundary of the region in which strange matter at the surface 
(p = 0) is absolutely stable (B l / A < 165 MeV). From this figure, we see that at a 
critical value of B 1 ^ = 165 MeV, B eS can be as large as 200 MeV inside the star for 
A > 4. Since such large values of B c r effectively soften the EOS of strange matter, 
the bag model with a pressure dependent bag constant allows for the possibility 
of strange stars that are more compact than those predicted by the bag model, as 
plotted in Fig. 6, for fixed B. 

In order to see how the resultant softening of the strange matter EOS affects 
the upper limit of the mass, M up , allowed by the inferred radiation radius of the soft 
X-ray source RX J1856. 5-3754, we plot the contours of M up in the B 1 / 4 -A plane (see 
Fig. 10). When A > 3, quark stars which are more massive than those allowed by the 
bag model prediction (A = 3) exist in the region of the absolute stability of strange 
matter at zero pressure. Because the inferred radiation radius of RX J1856. 5-3754 
is very small, however, the change is small for A ~ 4. 
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Fig. 10. Contours of the upper limit of the mass, M up , allowed by the inferred radiation radius 
of the soft X-ray source RX J1856.5-3754 in the B 1/4 -A plane, calculated for a c = 0. The 
contours of M up = IMq (solid curves), M up = 0.9Mq (dashed curves), M up — 0.7Mq (dotted 
curves), M up — 0.5Mq (long dashed curves) and M up = 0.3Mq (dot-short-dashed curves) are 
included. The vertical dot-dashed curve is the boundary of the region in which strange matter 
is absolutely stable at zero pressure (B 1 ^ < 165 MeV). 
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